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Abstract 

In this paper we study the class of backward doubly stochastic differential equations 
(BDSDEs, for short) whose terminal value depends on the history of forward diffusion. 
We first establish a probabilistic representation for the spatial gradient of the stochastic 
viscosity solution to a quasilinear parabolic SPDE in the spirit of the Feynman-Kac for- 
mula, without using the derivatives of the coefficients of the corresponding BDSDE. Then 
such a representation leads to a closed-form representation of the martingale integrand of 
BDSDE, under only standard Lipschitz condition on the coefficients. 
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1 Introduction 



Backward stochastic differential equations (BSDEs, for short) were firstly been considered in 
it linear form by Bismut [H [2] in the context of optimal stochastic control. However, nonlinear 
BSDEs and their theory have been introduced by Pardoux and Peng [12] ■ It has been enjoying 
a great interest in the last ten year because of its connection with applied fields. We can cite 
stochastic control and stochastic games (see [8]) and mathematical finance (see [6]). BSDEs 
also provide a probabilistic interpretation for solutions to elliptic or parabolic nonlinear partial 
differential equations generalizing the classical Feynman-Kac formula [EH E]. A new class 
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of BSDEs, called backward doubly stochastic differential equations (BDSDEs, in short), was 
considered by Pardoux and Peng [15]. This new kind of BSDEs present two stochastic integrals 
driven by two independent Brownian motions B and W and is of the form 

Y 8 = / f(r,Y r ,Z r )dr+ [ g(r, Y r , Z r ) j dB r 

J s J s 

-J T Z r dW r , se[t,T], (l.i) 

where £ is a square integrable variable. Let us remark that in the sens of Pardoux Peng, the 
integral driven by {B r } r > is a backward Ito integral and the other one driven by {W r }r>o 
is the standard forward Ito integral. Further, backward doubly SDEs seem to be suitable 
giving a probabilistic representation for a system of parabolic stochastic partial differential 
equations (SPDEs, in short). We refer to Pardoux and Peng pj)] for the link between SPDEs 
and BDSDEs in the particular case where solutions of SPDEs are regular. The general situation 
is much delicate to treat because of difficulties of extending the notion of viscosity solutions to 
SPDEs. The stochastic viscosity solution for semi-linear SPDEs was introduced for the first 
time in Lions and Souganidis [10]. They used the so-called "stochastic characteristic" to remove 
the stochastic integrals from an SPDE. Another way of defining a stochastic viscosity solution 
of SPDE is via an appeal to the Doss-Sussman transformation. Buckdahn and Ma 0, H] were 
the first to use this approach in order to connect the stochastic viscosity solutions of SPDEs 
with BDSDEs. 

In this paper we consider the approach of defining stochastic viscosity solution of SPDEs 
given by Buckdahn and Ma [H H] which, in our mind is natural and coincide (if g = 0) with 
the well-know viscosity solution of PDEs introduced by Crandall et al [5]. In this fact, we will 
work in the sequel of this paper with the version of backward doubly SDEs introduced in [HI [4], 
which is in fact a time reversal of that considered by Pardoux and Peng [15]. Indeed, for I, f 
be Lipschitz continuous functions in their spatial variables and g G C°' 2 ' 3 ([0, T] x IR d x IR; IR d ), 
they consider a class of backward doubly SDEs is of this following form: 

Y^ = l(Xl> x ) + f f(r, X?, Y^ x , Zjf) dr + f g(r, X^, Y?) dB r 
Jo Jo 

- [ Z\: x I dW r , s G [0,t]. (1.2) 
Jo 

The diffusion process X t,x is the unique solution of the forward SDE 

Xl' x = x + J 6(r, Xl' x ) dr + a(r, X?) [ dW r s G [0, t], (1.3) 

where b and a are some measurable functions. Here the superscript (t, x) indicates the depen- 
dence of the solution on the initial date (t, x), and it will be omitted when the context is clear. 
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Buckdahn and Ma proved in their two works [HI 0], among other things, that u(t,x) = Y t t,x is 
a stochastic viscosity solution of nonlinear parabolic SPDE: 

du(t, x) = [Cu(t, x) + f(t, x, u(t, x), (Vua)(t, x))] dt 

+g(t,x,u(t,x)) dB t , (t, x) G (0, T) x IR d , (1.4) 
u(0,x) = l(x), x G IR d , 

where C defined by 

^ n k n 

i,j=l 1=1 j=l 

is the infinitesimal operator generated by the diffusion process X t,x . More precisely, they show 
thank to the Blumenthal 0-1 law that 

u(t, x) = IE |/(Xf) + J f{r, X f r ' x , Y^, Z l r ' x ) ds 

+ ^g(r,X t r >*X' X )dB r \F t B \. (1.5) 

It is well know that u is a jF B -measurable field. However, to the best of our knowledge, to date 
there has been no discussion in the literature concerning the path regularity of the process Z 
when / and I are only Lipschitz continuous, even in the special cases where the coefficients are 
enough regular. 

Our goal in this paper is twofold. First we show that if the coefficients I and / are con- 
tinuously differentiable, then the viscosity solution u of the SPDE (jl.4p will have a continuous 
spatial gradient d x u and, more important, the following probabilistic representation holds: 

d x u{t,x) = E|/(X^)iV* + ^ S /(r,Xf X'^Kn^dr 

+ j\(r,X^X' X )K dB r | } (1.6) 

where iV s is some process defined on [0, s], depending only on the solutions of the forward SDE 
(11 .3p and its variational equation respectively. This representation can be thought of as a new 
type of nonlinear Feynman-Kac formula for the derivative of u, which does not seem to exist in 
the literature. The main significance of the formula, however, lies in that it does not depend on 
the derivatives of the coefficients of the backward doubly BSDE (jl.2p . Because of this special 
feature, we can then derive a representation 

= dtufaXyMttX**), s e [0,t], (1.7) 
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under only a Lipschitz condition on I and /. This latter representation then enables us to prove 
the path regularity of the process Z, the second goal of this paper, even in the case where 
the terminal value of Y t)X is of the form l(X to , ...,X tn ), where 7r : = t < .. < t n = t is any 
partition of [0, t], a result that does not seem to be amendable by any existing method. 

Let us recall that this two representations have already be given by Ma and Zhang [TT] in 
the case of a probabilistic representation for solutions of PDEs via BSDEs. Consequently our 
approach is be inspired by their works. However, there are particularities: first, the derivative 
notion is take in the flow sense (independent of uji) because u, the stochastic viscosity solution 
of the SPDE (II .41) . is a random field. Secondly, the proof to the continuity of the representation 
of the process Z need, since F* = (Ff ® ^Yt)o<s<t is not a filtration, G* = (Ff ® ^F^)o<s<t 
which is a filtration. 

The rest of this paper is organized as follows. In section 2 we give all the necessary pre- 
liminaries. In section 3 we establish the new Feynman-Kac formula between coupled forward 
backward doubly SDE (iTTTjl - lOl) and the SPDE (Ol) . under the C 1 - assumption of the coef- 
ficients. The section 4 is devoted to give the main representation theorem assuming only the 
Lipschitz condition of the coefficients I and /. In section 5 we study the path regularity of the 
process Z. 

2 Preliminaries 

Let T > a fixed time horizon. Throughout this paper {W t , < t < T} and {B t , < 
t < T} will denote two independent (i-dimensional Brownian motions defined on the complete 
probability spaces (Qi,^, Pi) and (f2 2 , IP2) respectively. For any process {U s , < s < T} 
defined on jFj, IPj) {i = 1, 2), we write T^ t = a(U r — U s , s < r <t) and = T§ v Unless 
otherwise specified we consider 

Q = fti x Q 2 , T = T X ®T 2 and P = Pi ® P 2 . 

In addition, we put for each t G [0,T], 

F = {JF S = Tf ® F™ T v M \0 < s <T} 

where Af is the collection of P-null sets. In other words, the collection F is P-complete but 
is neither increasing nor decreasing so that, it is not a filtration. Let us tell also that random 
variables u)\ G fl± and ((u>2), ^2 G fi 2 are considered as random variables on fl via the 

following identification: 

u 2 ) = £M; C(uj u uj 2 ) = CM- 

Let E denote a generic Euclidean space; and regardless of its dimension we denote (; ) to be 
the inner product and |.| the norm in E. if an other Euclidean spaces are needed, we shall label 



4 



them as E\\ E 2 , ., ., etc. Furthermore, we use the notation d t — J^, d x — ••, gf^) and 

d 2 = = {d 2 x . X] )l j=1 , for G [0,T] x ]R d . Note that if ^ = (V>\ : IR d -> IR d , then 

— (9xjtp l )fj =1 is a matrix. The meaning of (9^, (9 ra , etc. should be clear from the context. 

The following spaces will be used frequently in the sequel (let X denote a generic Banach 
space): 

1. For t G [0, T], L°([0, t\; X) is the space of all measurable functions ip : [0, t] i— > A\ 

2. For < £ < T, C([0, £]; A") is the space of all continuous functions ip : [0, t] \— > Af; further, 
for any p > we denote |v|o't = SU P llv 9 ! 5 )!!^ when the context is clear. 

0<s<t 

3. For any k, n > 0, C fc,n ([0, T] xi?; Ei) is the space of all Ei-valued functions </?(£, e), (t, e) G 
[0, T] x E, such that they are /c-times continuously differentiable in t and n-times contin- 
uously differentiable in e. 

4. C£([0,!T] x E; Ei) is the space of those ip G C^O,? 1 ] x £; £i) such that all the partial 
derivatives are uniformly bounded. 

5. Vl^ 1 ' 00 (£', £i) is the space of all measurable functions ip : E t— > such that for some 
constant > it holds that |t/>(x) — < -^1^ — y\E,^ x ,U £ -E 1 - 

6. For any sub-a-field £ C JF|? and < p < oo, £ p (<?; £") denote all E'-valued ^-measurable 
random variable £ such that IE|£| P < oo. Moreover, £ G L°°(Q;E) means it is im- 
measurable and bounded. 

7. For < p < oo, L P (F, [0, T}; X) is the space of all Af-valued, F-adapted processes £ 
satisfying IE 11611*^ < °°; and also, £ G L°°(F, [0, T]; IR d ) means that the process 
£ is uniformly essentially bounded in (t,u). 

8. C(F, [0, T] x i^) is the space of Ei-valued, continuous random field <p : Q x [0, T] x E 1 , 
such that for fixed e G E, </?(., ., e) is an F-adapted process. 

To simplify notation we often write C([0, T] x £7; Si) = C°'°([0, T] x £; E^; and if E x = IR, then 
we often suppress £ x for simplicity (e.g., C fe ' n ([0, T] x E; IR) = C fc ' n ([0, T]x£), C k ' n (F, [0, T] x 
-E; IR) = C k,n (F, [0,T] x E),..., etc.). Finally, unless otherwise specified (such as process Z 
mentioned in Section 1), all vectors in the paper will be regarded as column vectors. 

Throughout this paper we shall make use of the following standing assumptions: 

(Al) The functions a G C$'\[0,T\ x IR d ; M dxd ), b G C 6 O,1 ([0,T] x IR d ; IR d ); and all the partial 
derivatives of b and a (with respect to x) are uniformly bounded by a common constant 
K > 0. Further, there exists constant c > 0, such that 

fa(t,x)a(t,x) T Z > c|£| 2 ,Vx,£ G IR d , t G [0,T]. (2.1) 
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(A2) The function / G C{T B \ [0,T] x IR d x IR x IR d ) n W^^T] x IR d x IR x IR d ) and 
/ £ W 1,0 °(1R d ). Furthermore, we denote the Lipschitz constants of / and I by a common 
one K > as in (Al); and we assume that 

sup {|6(f,0)| + \a(t,0)\ + 0,0)| + |«?(0)|} < K. (2.2) 

0<t<T 

(A3) The function g G C°' 2,3 ([0, T] x IR d x IR; M d ) 

The following results are either standard or slight variations of the well-know results in 
SDE and backward doubly SDE literature; we give only the statement for ready reference. 

Lemma 2.1 Suppose that b G C(F, [0, T] x M d ; M d ) n L°(F, [0, T]; W 1 '°°(R d ; M d )), 
a G C(F, [0, T] x M 1 - m dxd )r)L°(F, [0, T];W 1 '°°(R d ] IR dxd )), with a common Lipschitz constant 
K>0. Suppose also that b(t,0) G L 2 (F, [0,T}; IR d ) and a(t, 0) G L 2 (F, [0, T]; IR dxd ). Let X be 
the unique solution of the following forward SDE 

X s = x + J b{r, X r ) dr + J a{r, X r ) dW r . (2.3) 
Then for any p>2, there exists a constant C > depending only on p,T and K, such that 

E(\X\$) < C | \x\* + Je[\\Ks, 0)| p + W(s, 0)H ds} (2.4) 

Lemma 2.2 Assume f G C(F, [0,T] x Mx M d )nL°{F, [0, T]; W 1 ' 00 (iR x i? d )) ; witfi a uniform 
Lipschitz constant K > 0, such that f(s, 0, 0) G L 2 (F, [0,T]) and 5 G C(F, [0,T] x JR x 
iR d ; 5?*) n L°(F, [0,T]; W 1 ' 00 (iR x i? d ; iR z )) with a common uniform Lipschitz constant K > 
respect the first variable and the Lipschitz constant < a < 1 which respect the second 
variable and such that g(s, 0,0) G L 2 (F, [0, T]). For any £, GL 2 (jF ;iR), let (Y, Z) be the adapted 
solution to the BDSDE: 

Y s = £+ [ f(r,Y r ,Z r )dr+ [ g(r,Y r , Z r ) dB r - [ Z r I dW r . (2.5) 
Jo Jo Jo 

Then there exists a constant C > depending only on T and on the Lipschitz constants K and 
a, such that 

e£ \Z s \ 2 ds < CIE^\ 2 + j\\f{s, 0, 0)| 2 + \g(s, 0, 0)| 2 ] tfe j . (2.6) 
Moreover, for all p > 2, there exists a constant C p > 0, such that 

IE(\y\*o P t) < CpIE^ + £[\f( 8 , 0, 0)1* + \g(s, 0, 0)H ds} (2.7) 
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We now review some basic facts of the anticipating stochastic calculus, especially those related 
to the backward doubly SDEs (see Pardoux-Peng [T5]). For any random variables £ of the form 



Z = F[ (pidWt,.., (f n dW s ; fadB.,..., ^j p dB s 
\Jo Jo Jo Jo 

with F G CZ°(M n+p ), <pi,...,<p n G L 2 ([0,T],IR d ), e L 2 ([0,T],IR d ), we let 

™ <9F / /" T /" T /" T \ 

D ^ = Y1q^\J VidW t ,..,J <PndW s ;j ^dB,, J i/; p dB a J <pi(t). 

For such a £, we define its 1, 2-norm as: 

||£|| 2 , 2 = IE |£| 2 + 1E / |L» r £| 2 dr 

5 denoting the set of random variable of the above form, we define the Sobolev space 

D l,2 A 

The "derivation operator" D, extends as an operator from ID 1 ' 2 into L 2 (f2, L 2 ([0, T], IR d )). 

We shall apply the previous anticipative calculus to the coupled forward backward doubly 
SDEs (jl.3P - (jl.2p . In this fact, let us consider the following variational equation that will play 
a important role in this paper: for % = 1, .., d, 



t,x 



e t + I d x b(r,Xjr B )V i X^ B dr + ^2 [ d^^X^ViX^ J, dW?, 

Js j=1 Js 



+ [\d x f(r, S*^(r))Vi^ + d y f(r, ~^(r))V^ + (d z f(r, S*^(r)), V<#">]dr 

+ r[^(r,0^(r))V J X^ + ^(r,e^(r))V^]rf J B r - f [ dW r , (2.8) 

Jo «/ 

where ^ = (0, 1, 0) T G H d ,H t ' x = (9^, Z*'*), = (X t>x , Y t>x ) and a j (.) is the j-th 
column of the matrix cx(.). We recall again that the superscription i)X indicates the dependence 
of the solution on the initial date (t, x), and will be omitted when the context is clear. We also 
remark that under the above assumptions, 

(VX t ' x ,VY t ' x ,VZ t ' x ) G L 2 (F;C([0,T];IR dxd ) x C([0,T];IR d ) x L 2 ([0, T]; M dxd )). 

Further the d x d-matrix- valued process VX*' Z satisfies a linear SDE and VX l t ,x = I, so that 
[VX*' X ] _1 exists for s G [0, t], IP-a.s. and we have the following: 
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Lemma 2.3 Assume that (Al) holds; and suppose that f G C b ' ([0, T] x M 2d+1 ) and 

g G C°' 2,3 ([0,T] x m d+l ;IR d ). Then (X,Y, Z) G L 2 ([0, T]; E) 1 ' (IR 2d+1 )) , and there exists a 

version of (D s X r , D s Y r , D s Z r ) that satisfies 

D s X r = VX r (VX s )- 1 ( x(s,X s )l {s < r} , 

D s Y r = Vr r (VX s )- 1 a(s,X s )l {s < r} , 0<s,r<t. (2.9) 
D s Z r = VZ r (VX s )-V(s,X s )l {s < r} , 

Lemma 2.4 Suppose that F G iD 1,2 . Tnen 

(i) (Integration by parts formula): for any u G Dom(S) such that Fu G L 2 ([0,T] x f2; one 

has Fu G Dom(5), and it holds that 

[ (Fu t , dW t ) = 5{Fu) =F f (u t , dW t ) - [ D t Fu t dt; 
Jo Jo Jo 

(ii) ( Clark-Hausman-Ocone formula): 

F = IE{F) + f IE{D t F | F t }dW t . 
Jo 



3 Relations to stochastic PDE revisited 

In this section we prove the relation (11. 71) between the forward backward doubly SDE fll.2p - fll.3l) 
and the quasi-linear SPDE (11.41) . under the condition that the coefficients are only continuously 
differentiable. Indeed, since Buckdahn and Ma [HI 0] provide that, if / and I are only Lipschitz 
continuous, the quantity u(t, x) = Y^' x is a stochastic viscosity solution of the quasi-linear 
SPDE (11.41) . relation in (11.71) becomes questionable. Our objective is to fill this gap in the 
literature and to extend the results of Ma and Zhang [H] given in the case of the probabilistic 
interpretation of PDEs via the BSDEs. 

Theorem 3.1 Assume (Al) and (A3) and suppose that f G C^\[0,T] x M d x Ut x M d ) and 
I G Cl{IR d ). Let (X t,:r , Y l ' x , Z l ' x ) be the adapted solution to the FBDSDE O^-fTOL and set 
u{t,x) = Y t ' x the stochastic viscosity of SPDE Then, 



(i) d x u(t,x) exists for all (t,x) G [0,T] x IR d ; and for each (t,x) and i=l,...,d, the following 
representation holds: 

d Xi u{t,x) = E{dJ(X^)ViX^ x 

+ [ [d x f(r, E^WViX? + d y f(r, E t ^{r))V l Y^ + d z f(r, E f > x (r))V \Z^\dr 
Jo 

+ jT [d x g(r, e^(r))V,X^ + d y g(r, Q^{r))Vi^]dB r | J (3.1) 
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where Q*' x = (X^',F^'), H*« a = (6*' x , Z*' x ), and (VX*> X , VF*'*, VZ*> X ) i/*e untgue sotorc of 
equation (\2.8i) : 

(ii) d x u(t,x) is continuous on [0, T] x M ; 

(Hi) Z l s ' x = d x u(s,Xl> x )a(s 7 Xl> x ),V s G [0,t], IP-a.s. 

Proof. For the simple presentation we take d — 1. The higher dimensional case can be treated 
in the same way without substantial difficulty. We use the simpler notations l x , (f x , f y , f z ), (g x , g y , g z ) 
respectively for the partial derivatives of I, f and g. 

The proof is inspired by the approach of Ma and Zhang |H] (see Theorem 3.1). Nevertheless, 
there exists slight difference due in the fact that the solution of SPDE's is a random field; more 
precisely will show that it is a conditional expectation with respect the filtration (J r ; B )o<t<T- 

We first prove (i). Let (t, x) G [0, T] x IR be fixed. For h ^ 0, we define: 

vt.x\h ~yt,x -yt,x+h ~yt,x yt,x+h yt,x 

vx h = ±s s_. V yh = _s s_. vz h = _^ s_ s e r Qj f ] 

h h h 

It follows analogously of the proof of Theorem 2.1 in pj)]) that 

TE{\AY%' 2 t = m{\VY h - VY l ' x \* j} -> as h -> 0. (3.2) 

We know also that processes Y t,x , Y t)X+h ^ VY h and AY h are all adapted to the a-algebra 
F* = (jF*) < s < t , where T\ = J-f ® J~7t- I n particular, since W is a Brownian motion on 
(ft 2 , J~2i IP2), applying the Blumenthal 0-1 law (see, e.g, |9J), Y^' x = u(t,x), Y^' x+h = u(t,x + 
h), VY t h = j^[u(t, x + h) — u(t, x)] and AY t h are all independent of (or a constant with respect 
to) u>2 G Q 2 . Therefore we conclude from the above that d x u exist, as the random field and 
d x u(t,x) = VYt' x , for all (t,x). Finally, taking the conditional expectation on the both sides 
of (12.80 at s = t, the representation (13 . 1 1) hold and finish the prove of (i). 

We now prove (ii). Let (U, Xi) G [0, T] X IR, i = 1, 2. Knowing that t\ and t 2 played inverse 
roles one another, we assume without losing a generality that t\ < t 2 . Since d x u is a conditional 
expectation with respect the filtration (JF^) < s < f , we have 

\d x u(t 1 ,x 1 )-d x u(t 2 ,x 2 )\ < IE{A(ti,xi) -A(t 2 ,x 2 ) I Fg} 

+ \TE{A(t 2 ,x 2 ) \rZ}-W,{A(t 2i x 2 ) (3.3) 

where 

A(t,x) = l x (X^ x )VX^ 

+ [\f x (r, *(r))VI^ + f y (r, E*'* (r))V^ + f z (r, ~<' x (r))VZf] dr 
Jo 

+ / [g x (r^ x (r))VXl' x + g y (r,Q^(r))VY^]dB r . (3.4) 
Jo 
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Thanks to the quasi-left-continuity of (J r ^)o< s <t, we see that 



lim \jE(A(t 2 ,x 2 ) | - m(A(t 2 ,x 2 ) \ = 0, (3.5) 

tllt2 

independently of x 2 . In virtue of (13.31) and (13.51) ). to prove {%%) it remain to show that 

lim JE{A(t 1 ,x 1 ) - A{t 2 , x 2 ) | Tf\ = 0. (3.6) 

To this end, since A(t, x) is a stochastic process and in virtue of Kolmogorov-Centsov Theorem 
(see [9]), it suffices to show that 

IE (\A(t 1} xx) - A(t 2 , x 2 )\ 2 ) < C(\h -t 2 \ 2 + \x x - x 2 \ 2 ), 

what we do now. Recalling the definition of A(U, Xi), % — 1, 2 and denoting 



and 



we get 



G*(r) = f x (r, S^(r))V^ + f y (r, &*(r))VY? + f z (r, S*-»(r))VZ, 



H^ x (r) = g x (r, &> x {r))V X l r ' x + g y (r, '(r))V^ 



t,x 



{AfaxJ - A(t 2 ,x 2 )\ < UX^VX^ 1 - l x (X**>)VX*>°»\ 



+ / \G tl,xl (r) -G t2 ' X2 (r)\dr + 



rti 

/ {H tx ' xi {r)~ H t2 > X2 (r))dB r 
Jo 



+ f 2 \G t2 ' X2 {r)\dr+ f 2 H t2 ' X2 {r)dB, 



'ti Jti 

Taking the expectation, it follows by Holder's and Burkolder-Gundy Davis inequalities that 

IE (\A{t 1>Xl ) - A(t 2 ,x 2 )\ 2 ) < CTEilUX^VX^ - l x {X l 2 ' X2 )V X l 2 ' X2 \ 2 



+ 



\G^ Xl (r) 



I' \H^ Xl {r) - H t2 ' X2 {r))dB r 
Jo 



G h ' X2 {r)\ 2 dr + 

+ f 2 \G t2 ' X2 (r)\ 2 dr + f 2 H t2 > X2 (r)dB r 
Jti Jti 

< cmimx^vxl 1 '* 1 - i x (x t2 ' X2 )wx t2 ' X2 \ 2 

+ f 1 \G tl ' Xl (r)-G t2 > X2 (r)\ 2 dr+ I ' \H tuXl (r) - H t2 ' X2 (r)\ 2 dr 
Jo Jo 

+ {t 2 -t 1 ) ! 2 IG^ir^dr+ih-h) [ * \H t2 > X2 (r)\ 2 dr\ . 

Jti Jtj J 

By similar standard computations in Ma and Zhang [11] (see proof of Theorem 3.1), we obtain 

IE (\A{t u x x ) - A(t 2 ,x 2 )\ 2 ) < C(\t 2 - t x \ 2 + \x 2 - Xl \ 2 ) 
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that provide the proof of (ii). 

It remains to prove (Hi). For a continuous function <p, let us consider {<f e } e >o a family of C 0,o ° 
functions that converges to ip uniformly. Since b, cr, I, / are all uniformly Lipschitz continuous, 
we may assume that the first order partial derivatives of b £ , a £ , I s , f £ are all uniformly bounded, 
by the corresponding Lipschitz constants of 6, a, I, f uniformly in e > 0. Now we consider the 
family of FBDSDEs parameterized by e > 0: 

x t,x = x + f b e (r, Xl' x )dr + /* a £ (r, X**) j dW r ; 

(3.7) 

= l £ (X^) + j; f £ (r, X?, Y?°, Z^)dr + j; g(r, X?, Y^)dB r - $° \ dW r 

and denote it solution by (X t ' x (e), Y t ' x (e), Z t,x (e)). We define u e (t,x) = Y^ x (e). Theorem 3.2 
of [15] provide that u £ is the classical solution of stochastic PDE 

du £ (t, x) = [C £ u(t, x) + f £ (t, x, u £ (t, x), (Vu £ a £ )(t, x))] dt 

+g(t, x, u £ (t, x)) dB t , (t, x) e (0, T) x M d , (3.8) 

u £ (0,x) = l £ (x), x e IR d . 

For any {x £ } C IR n such that x £ — > x as e — > 0, define 

(X e , y e , Z e ) = (X t,x " (e), Y t,xB (e), Z t,x *' '(e)). Then it is well know according the work of Pardoux 
and Peng pj)] that 



Y £ = u £ (s, X £ ); Z £ = d x u £ (s,X £ s )a £ (s,X £ s ), V s 6 [0, t], IP - a.s. (3.9) 
Now by Lemma 2.1 and Lemma 2.2, for all p > 2 it hold that 

IE - X|SJ + |F £ - Y\% + jf |Z S £ - Zfds} -> (3.10) 

as £ — > 0. Moreover let us recall (VX e , VF £ , VZ e ) the unique solution of the variational 
equation of (13.71) . Using again Lemma 2.1 and Lemma 2.2 we get 

e||VX £ - VX\* > p t + \VY £ -VY\* > p t + \VZ £ s - VZ s \ 2 ds^ 0, (3.11) 

as e — > 0. Thus it is readily seen that 

W>{l £ x (X £ )VX £ \F*} -> E{Z x (X )VX |jf }, 
F — > 0. Furthermore, by the analogue step used in [11], one can show that 



IE U [fl(r)VX £ + f y (r)VY £ + f £ (r)V 'Z £ }dr + j\g x (r)VX £ + g y (r)VY £ ]dB r \F t 
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converge to 



E | J [f x (r)VX r + f y {r)VY r + f z {r)VZ r ]dr + jf fe(r)VI r + ^(r)VF r ]^ r |^ f 

F £ — ► 0. Therefore, we get 

d x u £ (t,x £ ) — > d x u(t,x), as e — > IP — a.s., 
for each fixed (t, x) G [0,T] x IR. Consequently, possibly along a subsequence, we obtain 
Z £ = \imdu £ (s,X £ s )a £ (s,X £ ) = du(s, X s )a(s, X 9 ), ds x dP - a.e. 

e— »0 

Since for IP — a.e. u, d x u(., .) and X are both continuous, the above equalities actually holds 
for all s G [0,t], IP-a.s., proving {Hi) and end the proof. ■ 

The following corollary is the direct consequence of the Theorem I3.ll The convention on 
the generic constant C > still true. 

Corollary 3.2 Assume that the same conditions as in Theorem \3.1\ hold. and let (X t,x , Y t,x , Z l ' x ) 
be the solution of FBDSDE <\l.S^ - <\1.3il . Then, there exists a constant C > depending only on 
K, T, and for any p > 1, a positive (JF*) < S < 4 , IP)-process T l ' x , such that 

\d x u(t,x)\ < CT*'*, V (t, x) G [0, T] x IR d , IP-a.s. (3.12) 

Consequently, one has 

\Zl' x \ <CTf(l + |X*' x |), Vse[0,t], iP-a.s. (3.13) 

Furthermore, Vp > 1, t/zere exists a constant C p > 0, depending on K, T, and p such that 

n3{\X*\%+\Y*\%+\Z t *\%} < C p (l + \x\n. (3.14) 

Proof. We assume first that p > 2. The case 1 < p < 2 then follows easily from Holder 
inequality. By Lemma 2.1 and Lemma 2.2, we can find constant C > such that 

mS^VX t 'X' p t + \VY t '%' p t + ^ | VZl' x \ 2 dr^j \<C. 

Then, from the identity (13. II) . we deduce immediately that \d x u{t, x)\ < CY\ X , for all (t, a;) G 
[0,T] x IR, where 

rf = IE (\VX^\ + jf [|VJ^'*| + IVY^I + {VZ^Wdr + [VX l ; x + VY^ x ]dB r \ J^j . 
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Moreover we get for s G [0, t], IE(|r* ,a; | p ) < C. Then Theorem 13.11 (Hi) implies that 

\Zl' x \ < CTf (1 + \Xl' x \), Ws E [0, *], P-a.s. 

Now, applying again Lemma 2.1 and 2.2 and recalling (13. 13|) we get (I3.14p . for p > 2. ■ 

To conclude this section, we would like to point out that in Theorem 3.1, the functions 
/ and / are assumed to be continuously differentiable in all spatial variables with uniformly 
bounded partial derivatives, which is much stronger than standing assumption (A2). The 
following theorem reduces the requirement on / and I to only uniformly Lipschitz continuous, 
which will be important in our future discussion. 



Theorem 3.3 Assume (A1)-(A4), and let (X,Y,Z) be the solution to the FBDSDE fTOj) - 
(ji.ffi) . Then for all p > 0, there exists a constant C p > such that 

mi\X\** t + + ess sup \Z S A < C p (l + \x\ p ). (3.15) 

I 0<s<t J 



Proof. In the light of the corollary 13.21 we need only consider p > 2. By Lemma 2.1 and 
Lemma 2.2 it follows that for any p > there exists C p > such that 

m\X\% + \Y\$}<C P 0. + \x\*). (3.16) 

Next, by similar argument of Theorem 3.1 (Hi), we consider two sequences of smooth functions 
{f e } £ and {l e } e with their first order derivatives in (x, y, z) uniformly bounded in t and e such 
that 

lim<{ sup \f £ (t,x,y,z) - f(t,x,y, z)\ + sup \l £ (x) - l(x) 



(t,x,y,z) 



Denoting (X e , Y e , Z e ) the unique solution of the corresponding FBDSDEs and applying Corol- 
lary [321 we can nn d for any p > 2 a constant C p > 0, independent of e, such that 

E (TO) <C p (l + \x\v). (3.17) 

Furthermore, by (13.101) we know that IE / \Z\ — Z s \ 2 ds — > as s — > 0. Thus, possibly along 

Jo 

a sequence say (e n )n>i we have lim^oo Z £n = Z ds x cOP-a.s. Applying Fatou's lemma and 
recalling (\3.17h we the obtain 



E<jess sup \Z S \ P \ <C p (l + \x\ p < 

0<s<t 



which leads to (13.151) . as desired. 
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4 Representation theorem 



In this section we shall prove the first main theorem of the paper. This theorem can be 
regarded as an extension of the nonlinear Feynman-Kac formula obtained by Pardoux-Peng 
[T5] . It gives a probabilistic representation of the gradient (rather than the solution itself) 
of the stochastic viscosity solution, whenever it exists, to a quasi-linear parabolic stochastic 
PDE. Unlike the cases studied in d 3 . 1 H . in this section, our representation does not depend 
on the partial derivatives of the functions /, I and g. In this context such representation is 
the best tool for us to study the path regularity of the process Z in the BDSDE with non- 
smooth coefficients. For notational simplicity, we shall drop the superscript t,x from the solution 
(X, Y, Z) of FBDSDE (O)l -(j01). 

To begin with, let us introduce the two important stochastic integrals that will play a key 
role in the representation: 



( [a-\r,X T )VX T } T idW T 

J r 



and 



N° = — *— (M;) T [VX r ]-\ < r < s < t. 
s — r 

Let us recall that 

e|m;| 2p < CpmH \a-\T,x T )vx T \ 2 dT^ (4.i) 
< c p ( s - rfm (i vx T |:f ) < c p ( s - r y, 

where C p > is a generic constant. 

An other hand, let us define the filtration G* = {j 7 ^ ® ■, < s < t} which will play a 
important role in the proof of the continuity of the process Z in the BDSDE. 

Lemma 4.1 For any fixed t G [0, T] and any H G L 00 ^ 1 , [0, T]; M) we have 

(i) IE\ J^H r M r s dB r \ <+oo 



ii) for IP.a.e., u G fi, the mapping s i— > f* ■^H r {uj)Ml{uj)dB r {uj) is continuous on [0, t] 
Hi) for IP.a.e. , u G Q, the mapping s i— > IE{f* ■^H r MldB r /Q t s \{u) is continuous on [0,t] 



Proof. First, for any < r < s < t we denote 

/ T * ^- r H r M s r dr } 0<r<s 
AZ { (4.2) 
0, if s = t. 
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To simplify notation, when r = we denote A$ = A s . Further, let (3 be such that a = 1—2(3 < \ 
and (3 < 1. Consider the random variable 

K 2 I 



M* 



sup 



0<t!<t 2 <t (^2 



(4.3) 



then by fl4~2|) and Theorem 2.1 of Revuz-Yor [E], we see that IE[M*] < +oo. 



To prove (i) we note that for any 0<r<s<tby Burkholder-Gundy- Davis's inequality 
one has 



JE\At\ < CJE 



< CJE 



H T M s r 


2 A 


dr I 


s — r 





1/2 



|M S | 2 
1 rl -dr 



r 



\2a 



1/2 



< CJE 



< CJE 





2 \ 


dr J 







1/2 



dr 



- r) 2 P 



1/2 



M* 



C( s _ r) (i/2)-/3 ]E( || i/ || ooM * )? (4 . 4) 



where ||.|| denotes the norm of L°°([0, T]). Again letting C > be a generic constant depending 
only on (3 and T, we have 

JE\A S T \ < C{E|| J ff||^ } 1/2 {lE(M*) 2 } 1 / 2 

< C\\H\\ Lao{[0jT]xn) \\M*\\ L 2 in) < oo. (4.5) 

Setting r = in (14.51) we proved (i). 

To prove (ii) let r = and observe that, in view of (i), A s is a stochastic integral for 
< s < t. Consequently, the mapping s \— > A s is continuous on [0,t]. It remain to prove (Hi). 
In this fact, we remark that the right-hand side of the inequality (14.41) (with r = 0) is clearly 
in L 1 ; thus we check easily that the process A is uniformly integrable. Therefore, by similar 
step in Ma and Zhang [TTj (see proof for (Hi) of Theorem 4.1) it follows that the G*-optional 
projection of A, denoting °A S = JE(A s \Q t s ), s G [0, t], has continuous path. This prove (Hi), 
whence the lemma. ■ 

Theorem 4.2 Assume that the assumptions (A1)-(A4) hold, and let (X,Y,Z) be the adapted 
solution to FBDSDE (TOft-ATl). Then 

(i) the following identity holds IP-almost surely: 

Z s = E^l(X )N° + jf f(r,X r ,Y r ,Z r )N?dr + jf g(r,X r ,Y r )N?dB r \J^ a(s,X s ) 
V0<s<t; (4.6) 
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(ii) There exists a version of Z such that for TP-a.e. u G f2, the mapping s i— > Z a (u) is contin- 
uous; 

(Hi) If in addition the functions f and I satisfy assumptions of Theorem 3.1, then for all 
(t,x) G [0,T] x IR d it holds that 

d x u(t,x) = m{l(X Q )N t + J* f(s,X r ,Y r ,Z r )N t r dr + jf g(r,X r ,Y r )N 8 r dB r \T^ . 

(4.7) 

Proof. Again we shall consider only the case d — 1. We assume first that / G C^M) and 
/ G C^QO, T] x IR 3 ). Using the nonlinear Feynman-Kac formula of Pardoux and Peng [15] we 
obtain that for < s < t, 

u(s,X s ) = Y s = e|z(X„) + J f(r,X r ,Y r ,Z r )dr + J g(r,X r ,Y r )dB r \^y (4.8) 

Similar arguments to those used in the work of Ma and Zhang [H], provide the following: 

d x u(s,X s ) = JE^(X )N^ + J f(r,X r ,Y r ,Z r )N s r dr + J g(r, X r , Y r )N*dB r \J^ . 

In particular, setting s = t we obtain (14. Th . this proves (in). 

We now consider the general case. First we fix s G [0, £]. For <p — I, f, let ip £ G C°°,e > 0, 
be the mollifiers of ip, and let (Y e , Z e ) be the solution of the BDSDE in ( 11.21) with coefficients 
(I s , f £ ,g). Then for each e > 0, as the previous we get 

Zl = ]E^l%X )N^ + f%r,X,Y^,Z £ r )N:dr + g(r,X,Y r £ )N:dB r \^a(s,X s ). 

(4.9) 

Passing to limit as e goes to zero in (14.91) . we get (14.61) IP-a.s., for each fixed s G [0,t]. 

We should note that to prove part (i) we still need to show that (14 .6p actually holds for all 
s G [0, T], IP-a.s., but it is easy to see that this will follow from part (ii); that is, the process 
Z has a continuous version. Thus it remain to prove only (ii). To do this we first note that 

Z s = m^l(X )N* + J f(r,X r ,Y r ,Z r )N?dr + jf g(r ) X r ,Y r )N s r dB T \g t \a(s,X s ). 

(4.10) 

Lemma 4.1 in |11| and Lemma 4.1 imply that the mapping 



IE f(r, X r , Y r , Z r )N s r dr + J' g(r, X r , Y r )N s r dB r \Q, 
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is a.s. continuous on [0, t\. By the similar ideas used in Ma and Zhang [TT] replacing F* 
which here is not a filtration by G*, it follows that the mapping s i— > IE [1(X )Nq |C/*} is also 
continuous on [0,t]. Consequently, the right side of (I4.6P is a.s. continuous on [0, t], and hence 
(14. 6p holds for all s G [0,i], IP-a.s., proving (n), whence the theorem. ■ 



Remark 4.3 A direct consequence of Theorem 4-2 that might be useful in application is the 
following improvement of Theorem 3.3: assume that (Al) and (A2) hold, then for all p > 0, 
there exists a constant C p > depending only on T, K and p such that 



\x\ 



E{\X\$ + \Y\$+\Z\%} <C P (1 + 
Indeed, since by Theorem 4-1, Z has a continuous version, thus ( BZZS|) becomes (IX 1 1\ 



(4.11) 



5 Discrete function case 

Let us recall that we have proved in Theorem 4.2 that the process Z in the solution to the 
FBDSDE (11.3p -( iL~2l) has continuous paths, under the condition that the coefficients / and I are 
only uniformly Lipschitz continuous. While such a result is already an improvement of that of 
Pardoux and Peng [l5] . it still within the paradigm of the standard FBDSDE in the literature, 
to wit, the terminal condition of the BDSDE is of the form 1(Xq) (see also [IH]). In this section 
we consider the class of BDSDEs whose terminal conditions are path dependent. More precisely, 
we assume that the terminal condition of the BDSDE is the form £ = l(X to , X tl , X tn ), where 
= to < h < ■■■■ <t n = tis any partition of [0, t]. We shall prove a new representation theorem 
for the process Z, and will extend the path regularity result to such 

Theorem 5.1 Assume that (A1)-(A3) hold; and in (A3), I : R^ n+1 ) -> M, Let tc : = t < 

t\ < < t n = t be a given partition of [0, t], and let (X, Y, Z) be the unique adapted solution 

to the following FBDSDE: 

X s = x+ j b(r,X r )dr+ f a(r,X r )dW r , 

J s J s 

Y s = l(X t0 ,X h ,...,X tn )+ I* f(r,X r ,Y r ,Z r )dr (5.1) 

Jo 

ps ps 

+ g(r, X r ,Y r )dB r - Z r dW r , s e [0, t] . 
Jo Jo 

Then on each interval (tj_i,tj),i = 1, ....,n, the following identity holds: 

Z s = M{l(X t0 ,X tl ,...,X t jNl_ i + J* f(r,X r ,Y r ,Z r )N? vti _Jr 

+ jT g(r,X r ,Y r )N° VU i dB r \F^o-(s,X s ). s G (t^U) (5.2) 
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Further, there exists a version of process Z that enjoys the following properties: 



(i) the mapping s i— > Z s is a.s. continuous on each interval (tj_i,ij), i — 1, , n; 

(ii) limits Z t - = lim s + t . Z s and Z f + = lim s i t . Z s exist; 

i i 

(Hi) \/p > 0, there exists a constant C p > depending only on T, K and p such that 

E\AZ U | p < C p (l + \x\ p ) < oo. (5.3) 

Consequently, the process Z has both cadlag and caglad version with discontinuities to,...,t n 
and jump sizes satisfying (15.31) 

Proof. As before we will consider only the case d — 1, and we assume first that f, I G C£. 
Let us first establish the identity (15.21) . We fix an arbitrary index i and consider the interval 
(U-ii ti)- By virtue of the Malliavin operator D, Theorem 2.4 and the uniqueness of the adapted 
solution to BDSDE, we obtain 

Z s = Y^d 3 W s X tj +[ [f x (r)D s X r + f y (r)D s Y r + fz(r)D s Z r ]dr 

PS PS 

+ / [g x (r)D s X r + g y (r)D s Y r + g z (r)D s Z r ]dB r - / D s Z r dW r 
Jo Jo 

S^djlVXt. + [\f x (r)V s X r + f y (r)V l Y r + f z (r)V l Z,]dr 
. j>i Jo 

+ jT [g x (r)VX r + g y (r)V % Y r + g z (r)V l Z r }dB r - V%dW r | (VX s )~ 1 a(s, X s ) 
= V i Y s (VX s y 1 a(s,X s ), t^KsKU. (5.4) 
Taking the conditional expectation IEj.lJ 7 *} on two sides of (15.41) we obtain 

Z s = E|X>ZVX t . + J o [/,(r)V s X r + / y (r)V;F r + / z (r)V%]rfr|^|(VX s )-V( S ,X s ). 

(5.5) 

The rest of the proof is similar to the BSDE case. It is clear now that to prove the theorem 
we need only prove properties (i)-(iii), which we will do. Note that (i) is obvious, in light of 
Theorem 4.2 and thanks to representation (15 .2p . Property (ii) is a slight variation of Lemma 
4.1 and Lemma 4.1 of Ma and Zhang [TT], with there being replaced by tj_i, for each i. 
Therefore we shall only check (Hi). To this end, we define AZ ti = Z ti+ — Z ti _. From (15.41) it 
not difficult to check that 

V- = v%[vx t j-V(t,,x ti ) Z ti+ = V* +1 YUVXj-V(t 4 ,X ti ). 
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Denoting a\ = -(V l+1 Y S - V l Y s ), i = 1, n, we have 

AZ k = (V t+1 Y S - V%)a(t i} X ti ) = -a^U, X u ). (5.6) 
Further, since (V*Y, V l Z) denotes the adapted solution of the following BDSDE 

V l Y T = y^d j lVX tj + [ [f x (r)VX r + f v (r)V l Y r + f z (r)VZ r }dr 

+ / \g x (r)VX r + g y (r)V i Y r + g s (r)V i Z r ]dB r - [ \7 l Z r dW r , re[W], 
Jo Jo 

if we denote (3i s = — (V l+1 Z s — V l Z s ), then we have 

al = ^V ti + [ S [f y (r)4 + f z (rMdr+ T[^(r)< + fc (r)^]dS r 
Jo Jo 

fcdWr, s G [0,t]. (5.7) 



So (a*,/? 1 ) is the adapted solution to the linear BDSDE (15.71) . It follows by Lemma 2.2 that 
Wp > there exists a C p > such that IE{|aJ | p } < C p . On the other hand the same estimate 
holds for a(s,X s ) because of assumption (Al) and Theorem 3.3; for [VX] -1 since it is solution 
of a appropriated SDE. It readily seen that (15.31) follows from (15.61) which prove (in). 
Finally, we note that when / and I are only Lipschitz, (15.21) still holds, modulo a standard 
approximation the same as that in Theorem 4.2. Thus properties (i) and (it) are obvious. To 
prove (Hi) we should observe that the standard approximation yield that AZ t £ . — > AZ ti a.s. So 
if (I5.3P holds for AZ t £ ., then letting e — > 0, (15.31) remains true for AZ t ., according the Fatou's 
lemma; that end the proof. ■ 
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